The stress field in granular systems: Loop forces and potential formulation 
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The transmission of stress through a marginally stable 
granular pile in two dimensions is exactly formulated in terms 
of a vector field of loop forces, and thence in terms of a single 
scalar potential. 

The loop force formulation leads to a local constitutive 
equation coupling the stress tensor to fluctuations in the lo- 
cal geometry. For a disordered pile of rough grains (even with 
simple orientational order) this means the stress tensor com- 
ponents are coupled in a frustrated manner, analogous to a 
spin glass. 

If the local geometry of the pile of rough grains has long 
range staggered order, frustration is avoided and a simple 
linear theory follows. Known exact lattice solutions for rough 
grains fall into this class. 

We show that a pile of smooth grains (lacking friction) can 
always be mapped into a pile of unfrustrated rough grains. 
Thus it appears that the problems of rough and smooth grains 
may be fundamentally distinct. 
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It has long been recognized in engineering practice 
that granular materials which lack cohesion cannot be re- 
garded entirely as solids, nor (until 'fluidised' or yielding) 
can they be regarded entirely as fluids . In this Letter 
we explore the scenario that there exists a marginal state 
of granular matter in between solid and liquid. Most 
particularly, we are concerned with how that intermedi- 
ate state transmits stress. 

The marginal state is readily characterised in terms of 
the mean coordination number z for mechanical equilib- 
rium. If external loads, such as gravity, are exerted on 
a pile of N grains [|j then the Nz/2 intergranular forces 
must be able to adjust so as to achieve Nd{d+ l)/2 con- 
straints balancing the force and torque on each grain in d 
dimensions. For perfectly rough grains, where every con- 
tact supports friction, this gives a critical coordination 
number Zc = d+1, whilst for ideally smooth asymmetric 
grains with no friction Zc = d{d +1) Q, 

For mean coordination numbers z < z^ the pile cannot 
be stable and under a general loading it must rearrange 
or consolidate. By contrast for z > z^ the intergranu- 
lar forces are underdetermined by the conditions of force 
balance alone and the deformation of individual grains 
together with their local constitutive equations becomes 
relevant. In the first case we have granular motion and 
the pile is fluid or yielding, whereas in the latter case in- 
tragranular deformation is crucial and we have a solid of 
constitution influenced by that of the individual grains. 
It is interesting to note that for ideally smooth grains. 



Zc matches the maximum mean coordination attainable 
by sequential packing (and in two dimensions it is the 
absolute maximum possible j^), so a (random) pile of 
smooth grains is probably never solid. 

The marginal case z = Zc \s special in that the grains 
do not need to move and the intergranular forces, and 
hence distribution of stress, are determined by condi- 
tions of force balance alone. This has been identified 
as a paradigm problem of theoretical granular mechanics 
' ' [H • The macroscopic analogues of balancing force 
and torque are that the stress tensor a obey 



F 



A{(t) ^ 0, 
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where A{a) — ^{a ~ a"^) is the antisymmetric part of 
a and F is an external force field. We shall focus here 
on a region transmitting stress but otherwise unloaded. 
This system of continuum equations is incomplete with- 
out d{d— l)/2 further constitutive relations. The central 
mystery we are seeking to unravel is the nature of these 
relations (one relation in two dimensions) when they also 
arise from force balance at the granular level. 

The key idea of this Letter, where we limit ourselves 
for simplicity to the two dimensional case, is to focus 
attention on local loops around contact points and corre- 
sponding loop forces. Each local void, which we will label 
by an index I, is enclosed by a loop of grains and contact 
points; each contact point separates two voids and so is 
party to two loops. Around each loop we take a loop 
force fi to circulate in the anticlockwise direction, that 
is +// is contributed to each intergranular force around 
the loop with (notionally) positive sign in the anticlock- 
wise direction, and vice versa for the reaction forces. The 
resulting total force across a contact point is then a dif- 
ference between two contributing loop forces. 

The loop forces have two key features, the first being 
that they parameterise the intergranular forces in a man- 
ner that automatically satisfies balance of force at each 
contact point and on each grain. In terms of them we 
can write the stress tensor on each grain as 
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where rig is the vector connecting (anticlockwise) the two 
contact points shared by loop I and grain g and Ag is the 
grain area defined in fig. 1 |P . The vectors rig = 
form an anticlockwise loop round the contact points of 
loop I and likewise rig = forms a clockwise loop 
round the contact points of grain g and play a central 
geometrical role in our discussion. 
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The second key feature of the loop forces is that, being 
defined on the loops which are half as numerous as the 
grains, they are a comparatively coarse-grained quantity. 
The hidden nature of the constitutive equation can now 
readily be appreciated, because the only remaining con- 
straints that the loop forces have to obey is balance of 
torque around each grain, that is 

J2''^a^fi = (3) 

which is one equation per grain and hence two equations 
per loop. 

In essence the two torque equations per loop give us 
two macroscopic conditions, A{d') = plus a constitu- 
tive relation. To achieve this separation explicitly we 
postulate a smooth interpolation of the loop forces to a 
function of continuous position /(r), with each loop and 
grain having a nominal centre Ri and Rg, respectively. 
This also introduces a natural continuum already at the 
grain scale without going into elaborate coarse-graining 
schemes. Then the force moment on grain g is, to a first 
order Taylor approximation, 

Sg = Cg- V/ (4) 

where 

Cg^'^rigRig and Rig = Ri - Rg. (5) 

The geometric tensors Cg characterise the grain local 
geometry and their symmetric and antisymmetric parts 
have quite distinctive properties. It is readily shown that 

Cg = AgR + Pg (6) 

where R ~ (-"^i^o) ^'^^ VL'D\t antisymetric tensor (also 
corresponding to ^-rotation in the plane) and Ag is the 
area associated with grain g shaded in figure 1. The 
areas Ag have the convenient property of tessalating the 
plane. The symmetric part, Pg, reflects fluctuation in 
the departure of the local geometry from isotropy, being 
given by 

^9 = H \^^l9Slg - tlgtlg) <J) 
I 

where the vectors sig and Ug are shown on figure 1. It is 
a crucial point that neighbouring grains share equal and 
opposite contributions ss* and —tt to their P tensors, and 
from this it follows that the average {P)g of P over all 
the grains in any convex region G goes to zero at least 

— 1/2 

as fast as Nq ' , where Nq is the number of grains in 
G. Note that the averaging of the Pg does not depend 
on any assumptions of isotropy. For periodic lattices (P) 
vanishes over one unit cell. 

To underpin the significance of Cg as the relevant ge- 
ometrical order parameter, consider a region of area A 



whose boundary is defined by the contact points of the 
outermost grains. Averaging over Sg in the region using 
eq. (|) gives 

^=\Y.^9-^f^{Cg-Vf). (8) 

g 

Carrying out the same summation using eq. (|2|) we ob- 
serve that the sum over loops inside the region cancels 
out and the only contribution comes from the boundary 
vectors, r^, between boundary contact points, 

boundary 

where fb is the external loading on the region. This con- 
tour sum can be converted, using Stokes theorem, into 

a=jj dsR-Vf = R-{Vf) . (10) 

Since (Pn) = we identify (Cg) ~ AR and, on comparing 
to eq. (^, we obtain 

{Cg ■ V/) = (Cg) ■ (V/) . (11) 

This result is particularly reassuring because it yields the 
effective response of a macroscopic granular region to a 
force field in terms of an effective macroscopic character- 
istic property (Cg). It is equivalent to results in contexts 
such as disordered dielectrics and continuum elasticity 
that are derived by combining a field equation with ei- 
ther a constitutive relation or an energy functional. 

We can now construct the constitutive equation in 
terms of the continuum stress tensor, which is given from 
eq. (^0|) as (7 = Vf where V is exactly the curl operator 
in two dimensions. The condition that the stress tensor 
on each grain be symmetric then leads, over and above 
the condition that the continuous a be symmetric, to the 
requirement 

Tr(g • <t) = 0, (12) 

where Q = R P ■ R?^ is a rotated verion of the symmetric 
part of the local Cg tensor. Relation (|l^) is the contin- 
uum constitutive equation that provides the missing link 
between the stress and the local geometry. It has the 
striking feature that the coefficients Q are spatially fluc- 
tuating quantities that locally add to zero, so that any 
simple attempt to identify a non-vanishing mean field 
value {Q) leads back to Eq. (Q) and therefore yields no 
new constitutive information. In this sense the problem 
is analogous to spin glasses, in our case it being & that is 
subject to spatial couplings of random sign ]Tot . 

Constitutive equations coupling stress components in 
a local linear relation, represented here by ([T^), have 
been discussed before under the presumption that a non- 
vanishing macroscopic version of Q exists. Examples are 
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the hypotheses of "fixed principal axes" |Tl| (correspond- 
ing to Q being traceless) and "oriented stress linearity" 
that were introduced to describe sand piles . Hav- 
ing shown that the simple grain average of Q vanishes, 
we are therefore particularly motivated to consider any 
case where a non-zero effective Q can be identified. 

Drawing inspiration from spin problems, we focus on 
the analogy of an antiferromagnct in which the local P 
tensors have alternation of sign in an unfrustrated pat- 
tern. The key idea is that for systems where every loop 
has an even number of edges, we can label all grains -I- 
or - such that each grain is surrounded only by opposite 
sign neighbours in an unfrustrated way. Then we can 
distinguish locally averaged values P+ and P_ restricted 
to the respective -t- and — grains, where P+ + P- =0. If 
the pile has staggered geometrical order such that these 
are non-zero, then this leads to the identification 

Qeff = RP+R^ (13) 

The simplest such case where we can see this work is a 
periodic lattice (fig. 2) with a generally anisotropic unit 
cell comprising of two grains, which we duly label + and 
— . Using eq. (^) for each grain separately, we write down 
the mean and deviatory cell stresses, for which we find 

a^(s+ + 5_) = (ar + P++ P_) • V/ (14) 

and 

5^{s+-S-^ = [5AR + P+-P^Yvf , (15) 

where A is the unit cell area and 5 A — ^4+ — A^ is 
the difference between the areas associated with the two 
grains. Using the periodic properties of the lattice one 
readily obtains 

P+ = -P- = \Y1 ^kSk - tktk (16) 

k 

where the vectors Sk and tk are labelled on figure 2. Then 
the condition that both 5 as well as a be symmetric then 
leads to the constitutive equation ( [l^ ) with Q replaced 
according to (|l3|). 

The idea of a staggered order parameter can be further 
exploited to formulate a mapping from perfectly smooth 
to perfectly rough systems. Given a pile of perfectly 
smooth grains, we note that its behaviour is identical 
to the same pile of rough grains provided at each contact 
point we insert a vanishingly small perfectly rough ball 
bearing. It is readily checked that starting from z — 6 
for the smooth grains leads to z = 3 for the resulting 
set of rough grains (3/4 of which are the new bearings 
each having z = 2), so we duly arrive at a marginally 
connected set of rough grains. Furthermore, we auto- 
matically have an alternating distinction between orig- 
inal grains (-I-) and bearings (— ). Explicit evaluation 
shows that for each small grain, P_ cx e as bearing ra- 
dius e ^ 0, so in this case the appropriate geometrical 



order parameter is p = Limit^^oP+/e. This order pa- 
rameter, and hence the corresponding Qoff, can have a 
non-zero average for an orientationally ordered assembly 
of smooth grains, as computed explicitly for a periodic 
array |^ . This contrast to the case of general rough grain 
piles, which explicitly require staggered orientational or- 
der for a non-zero QeS, suggests that piles of perfectly 
smooth grains may be inherently distinct from piles of 
rough ones. 

Finally we introduce an alternative formulation of the 
problem where the loop forces are defined in terms of 
a potential field. Since each loop is associated with a 
two-dimensional force the number of degrees of freedom 
that the potentials should provide is twice the number 
of loops, which (up to relative correction of due 
to the boundary) is exactly the number of grains in the 
system. Thus, we require that the potential field provide 
one degree of freedom per grain leading to scalar grain 
potentials, il^g. A natural way to precisely define the 
potentials is to choose that the forces are derived from 
them via 

fi-^Y.^'igi^s . (17) 

and the potentials are then governed by eq. (|^). 

We now assume that the ■0-field can be continued 
throughout the system just as we did previously for the 
force field. Then expanding around the loop centres, we 
obtain 

fi - ^Ci ■ VV . (18) 

The geometrical tensor Ci = ^igRig is the loop ana- 
logue of Cg and its antisymmetric part is likewise asso- 
ciated with an area Ai, where these areas tessalate the 
plane. The symmetric part of Ci is sensitive to the choice 
of grain centre positions: as there are two grains per loop 
these constitute four degrees of freedom per Ci tensor 
and the expectation is that these tensors can thereby be 
chosen to be purely antisymmetric. Using C; — AiR sim- 
plifies the expression of the local stress in terms of the 
potentials to 

Sg = CgifVgVli, (l9) 

and the continuum stress becomes 

<7 = VgVi%lj curl curlT/- . (20) 

It follows that the continuous version of -0 is exactly the 
Airy stress function |^ . 

To conclude, we have formulated a theory of stress 
transmission in granular matter at the marginal state be- 
tween solid and fluid in terms of loop forces, which in turn 
can be defined in terms of a scalar potential. Granular 
systems can possess a predominant antiferromagnetic or- 
der but with local frustration, much like in antiferromag- 
netic Ising systems. An antiferromagnetic order yields a 
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FIG. 1. The local geometry around a grain g. The vectors 
fig connect contact points clockwise around each grain g and 
anticlockwise around each void loop I, whilst the vectors Rig 
connect from grain centres to loop centres. The stress field 
is coupled through tensors Cg = Rigfig, whose antisym- 
metric part has magnitude Ag which is the area shaded in 
the figure. The symmetric parts of these tensors evaluate to 
Pg = 1 sigSig ~ UgUg KvA play a crucial role in the con- 
stitutive equation. As can be seen from the labelling of the 
vectors sjg and tig in the figure, neighbouring grains share 
equal and opposite contributions to their repsective Pg, so 
these tensors tend to zero upon local averaging. 



global consitutive relation leading to a macroscopic the- 
ory while frustration requires a more sophisticated ap- 
proach. The similarity between random granular piles 
and spin glasses suggests a possible use of methods from 
the latter to this problem. We have mapped smooth piles 
to antiferromagnetically ordered perfectly rough ones ob- 
taining for them an exact constitutive relation. This re- 
sult points to an inherent difference between rough (frus- 
trated) and smooth (ordered) systems and may thus un- 
dermine the idea of gaining insight into behaviour of piles 
of finite friction grains from understanding the perfectly 
smooth and perfectly rough limits. 

An extension of this work to three dimensions is cur- 
rently under way. 
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FIG. 2. Part of an anisotropic periodic lattice, where 
there are two distinct grains, labelled -I- and -, per 
unit cell. Stress transmission is controlled by the tensor 
P-i- — — P- = \^^y,8h8h — thtk- The lattice vectors are 
ii + S3, t2 + Si and ta -I- 5*2 and because these are not in- 
dependent the value of P+ turns out to be independent of the 
choice of loop centre. 



and kinematics of granular materials (Cambridge Univer- 
sity Press, Cambridge 1992); V. V. Sokolowskii, Statics 
of granular materials (Pergamon Press, Oxford 1965); 
A grain is defined as a rigid particle of finite volume and 
arbitrary shape. 

S. F. Edwards and D. V. Grinev, Phys. Rev. Lett. 82, 
5397 (1999). 

R. C. Ball, in Structures and dynamics of materials in 
the mesoscopic domain, eds. M. Lai, R. A. Mashelkar, 
B. D. Kulkarni and V. M. Naik (Imperial College Press, 
London 1999). 



[2] 
[3] 
[4] 

[5] 
[6] 

[7] 

[8] 



:ond-mat/9810124 



R. C. Ball and D. V. Grinev, 

RGB acknowledges correspondence with M. E. Gates on 
this point. 

S. F. Edwards and R. B. Oakshot, Physica D 38, 88 
(1989). 

See eg, S. F. Edwards and G. C. Mounfield, Physica A 
226, 1, 12, 25 (1996); J. -P. Bouchaud, M. E. Gates and 
P. J. Claudin, J. Phys. II (France) bf 5, 639 (1995). 
[9] The definition of Ag as a grain area is useful because the 

[10] 



[1] R. L. Brown and J. C. Richards, Principles of powder me- 
chanics (Pergamon, NY 1966); R. M. Nedderman, Statics 



entire system can be uniquely mapped in terms of these 
areas and Asystem ~ ^ -^g- 

The conclusion that the stress couplings have zero aver- 
age appears to be at variance with that of ^ who arrived 
at an equation similar (but not identical) to equation 
(p^. However, we find that their coupling tensors 
have similar averaging properties to our Pg. 
[11] J.P. Wittmer, P. Claudin, M.E. Gates and J.-P. 

Bouchaud, Nature 382, 336 (1996). 
[12] J.P. Wittmer, M.E. Gates and P. Claudin, J. Phys. I 
(France), 7, 39 (1997); M. E. Gates, J. P. Wittmer, J.- 
P. Bouchaud and P. Claudin, Phys. Rev. Lett. 81, 1841 
(1998). 

[13] N. I. Muskhelishvili, Some Basic Problems of the Math- 
ematical Theory of Elasticity, (P. Noordhofi', Groningen, 
Netherlands 1963). 



4 



